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Oh' Abstract. When subject to a pair-breaking perturbation, the pairing susceptibihty 

of a disordered superconductor exhibits substantial long-ranged mesoscopic fluctua- 
tions. Focusing on a thin film subject to a parallel magnetic field, it is proposed 

Cu ■ that the quantum phase transition to the bulk superconducting condensate may be 

preempted by the formation of a glass-like phase with multi-fractal correlations of a 
complex order parameter. Although not universal, we argue that such behavior may 

Pj I be a common feature of quantum critical phenomena in disordered environments. 
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PACS numbers: 74.55. H-h,74.62.-c,74.76.-w,73.21.-b 
> ■ Several years ago, in an inspiring sequence of papers, it was proposed by Spivak 

and Zhou jT| |2] that mesoscopic fluctuations can significantly influence the nature 
of the transition to superconductivity in the vicinity of the upper critical field, ifc2- 
(3 ■ Considering a magnetic field applied perpendicular to a film, it was argued that the 

^ ■ transition to the mixed phase is mediated by the formation of a superconducting 'glass- 

like' phase realized through the random Josephson coupling of droplets or domains 
nucleated at fields in excess of Hc2 [I]. Subsequently, focusing this time on the properties 
of a superconducting thin-film subject to a parallel field, qualitatively similar conclusions 
were drawn |2j. Recently, Galitski and Larkin [S] extended these ideas exploring the 



a. 

O ■ interplay of the proximity effect and quantum phase fluctuations on the integrity of the 

superconducting glass in the perpendicular field geometry. Here, focusing on a weakly 



% 



disordered thin-film subject to a parallel magnetic field H, we will offer a different 
perspective on the character of the superconducting phase close to the quantum critical 
point. In this geometry, taking the film thickness d to be smaller than the penetration 
depth, the field lines enter the sample and effect a mechanism of pair-breaking leading to 
the gradual suppression of the superconducting order parameter. However, in contrast 
to Ref. [2], we will find it convenient to limit considerations to the range d ^ gL^F^o/^, 
where gi denotes the Lande (^-factor, Xp is the Fermi wavelength, i the elastic mean free 
path, and ^o = yD/Ao (with the classical diffusion constant D = vpi/^) represents the 
diffusive superconducting coherence length of the unperturbed system. In this limit, 
where the impact of Zeeman splitting can be safely neglected, the transition to the 
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(gapless) superconducting phase is second order and described by a mean-field theory 
of Abrikosov-Gor'kov (AG) type [5]. 

Starting with a microscopic BCS Hamihonian for the disordered, symmetry-broken 
system, one can present the quantum partition function for an individual realization 
of the impurity potential as a functional field integral Z = f D[A,A] e~^^, where, to 
quadratic order in A, 



S^ = Y,j dvdv'A^^iv) 
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Here, Abcs denotes the coupling constant of the BCS interaction while nj^^(r,r') = 
^Ee„'^^n(r,r')G'c^m-e„(r,r'), where G,„(r,r') = (r|(«e„-i^)-V) is the sample specific 
Matsubara Green function, represents matrix elements of the (generically complex) 
Hermitian pairing susceptibility. Since we are primarily interested in the character 
of the transition, and not properties deep in the superconducting phase, the influence 
of the non-linear terms in the Hamiltonian can be neglected. 

Leaving aside the influence of dynamical fluctuations of the order parameter, in the 
saddle-point approximation, the transition to the superconducting phase is signaled by 
the appearance of solutions of the linear equation 

c?r'no(r,r')x(r')=z/7X(r), (2) 

with eigenvalues 7 > (Abcs^)^^) deflning z/ as the density of states (DoS) per spin 
of the normal phase (assumed constant). In the leading approximation, the influence 
of disorder on the symmetry-broken system can be explored by replacing the pairing 
susceptibility by its impurity average, (n)y. In the specifled thin-fllm geometry, the 
'paramagnetic term' due to the parallel magnetic fleld is negligible, and (II) y takes the 
form 

where 1/th = De'^{Hd)^/6 denotes the pair-breaking rate associated with the orbital 
motion of the particle in the external fleld |3]. In this approximation, when substituted 
into (j21), one flnds that the system becomes unstable against the formation of a 
homogeneous condensate when 1/th fulflUs the AG condition [HI, ln(Tc/T°) = ip{l/2) — 
-0(1/2 + 1/2'7ithTc), where T° is the transition temperature of the unperturbed system, 
and ip{z) = r'{z)/r{z) denotes the digamma function. In particular, when the pair- 
breaking rate equals to the order parameter of the unperturbed system, 2/ffj = Aq = 
2 X 0.88 T°, the bulk superconducting phase is destroyed altogether. 

However, as emphasized by Spivak and Zhou |2], mesoscopic fluctuations of the 
pairing susceptibility can signiflcantly influence the character of the transition. In 
particular, at the transition, the condensate wave function may acquire a texture in 
both phase as well as amplitude: Focusing on the static component, the complex 
pairing susceptibility exhibits spatial mesoscopic fluctuations 5Ho(r,r') = Ho(r, r') — 
(Ho(r,r'))v. For a given conflguration of the quenched impurity potential, a system 
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may lower its energy by forming a condensate where the corresponding wave function is 
random and moreover, because of the presence of the magnetic field, complex. In this 
case, an arrangement of supercurrent loops can optimally screen the applied field. In 
particular, the condensate can exploit the quenched spatial fluctuations of the pairing 
susceptibility to initiate a transition to a superconducting phase ahead of that predicted 
by the mean-field estimate above. However, in the bulk system, such an arrangement of 
supercurrent loops can be tolerated only in the near vicinity of the transition. Away from 
the quantum critical point, the energy cost in suppressing the amplitude of the order 
parameter in the superconducting phase becomes prohibitively high and a condensate 
of uniform phase must develop. 

To explore the influence of mesoscopic fluctuations on the nature of the transition 
one can proceed along complimentary paths: Firstly, one can develop a scheme 
perturbative in fluctuations 611, treating typical field configurations of the order 
parameter Ao(r) associated with the impurity averaged susceptibility as a variational 
Ansatz. When the scattering rate 1/th is in excess of 1/ffj, the superconducting order 
parameter Ao(r) exhibits spatial fluctuations around a zero mean with a correlation 
length Ln = ^oiffj/Tn — 1)^^^^ diverging at the critical point. Regarding these 
spatial configurations as domains or "droplets" of ordered phase with a characteristic 
length scale Ljy, fluctuations 6Ilo impose a long-ranged, complex, random Josephson 
coupling — a superconducting "gauge glass" |2|. While such a variational approach is 
naively applicable close to the bulk transition — see below — one might expect it to 
underestimate the impact of optimal fluctuations of the random susceptibility. 

Alternatively, treating the optimal fluctuations more accurately, one can pursue a 
more direct 'mean- field' approach, seeking explicit solutions of the stochastic saddle- 
point equation (0) (cf. Ref. [Sj). In this case, guided by the intuition afforded by 
the properties of band tail states of random Hamiltonian operators, one might expect 
the behavior close to the critical point to be dominated by spatially localized field 
configurations of Ao(r) associated with rare or optimal fluctuations of the random 
operator Hq. Ignoring the potential impact of spatial and dynamical fluctuations of 
the order parameter around these saddle-point configurations, such an approach would 
suggest that the formation of the homogeneous superconducting condensate phase is 
preempted by the nucleation of domains of ordered phase, the implication being that 
the bulk transition will be mediated by phase ordering of Josephson coupled islands. 
Crucially, since the droplets reflect fluctuations of the complex pairing susceptibility, 
the condensate wave function would itself be complex — a superconducting 'glass-like' 
phase of the type described by Ref. [2] . 

The distinction between the two approaches may to some extent be semantic. While 
the former approach may underestimate the capacity of the system to form 'droplets' of 
condensate, the nature of the bulk transition may rely more sensitively on the properties 
of the gauge glass: the geometry and rigidity of the superconducting domains is less 
important than the random Josephson coupling between the domains arising from the 
fluctuations. However, in the parallel field geometry, when taking into account the 
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long-ranged nature of the correlations of IIo! we will argue that the transition affords a 
potentially different interpretation. Specifically, we will show that fluctuations (5no(r, r') 
engage long-ranged diffusion modes and, as such, are correlated over long distances. Such 
long-range correlations discriminate the spectral properties of the pairing susceptibility 
from those of usual short-ranged random impurity models. In particular, in contrast to 
band tail states in low-dimensional Hamiltonian systems, the transition between the tail 
state region and the 'extended' bulk state region is not mediated by a region of weakly 
localised states but rather it is abrupt: Even states close to the band edge of IIo are 
power-law extended. If the condensate acquires the texture of these states, they will in 
turn inherit 'critical' spatial correlations. 

To explore in detail the spectral properties of the (static) pairing susceptibility 
Ho, one requires information about the full distribution function of matrix elements, 
P[no]. Clearly, being of two-particle character, such a programme seems unfeasible 
— at least in the two-dimensional geometry considered here. Instead, we will choose 
to characterize the scale of fluctuations through the variance of Hq which can be 
estimated diagrammatically. Here, to be consistent, one must discard fluctuations 
5no(r, r') occurring on a length scale |r — r'| shorter than the magnetic diffusion 
length Lh = \/Dth. Such fluctuations, which effect equally matrix elements of the 
unperturbed superconducting system, demand a more careful consideration of the 
influence of disorder on the pairing interaction itself. These effects are irrelevant 
to, and reach well-beyond, the present discussion. Although the ensemble average 
has all but vanished, on length scales |r — r'| ^ Lh-, one finds that the matrix 
elements of the pairing susceptibility become complex and exhibit significant long- 
ranged fluctuations. Taking into account the leading contribution from the impurity 
diagrams depicted in Fig. ^, a diagrammatic estimate of their magnitude obtains 
(5no(ri, r2)(5no(r3, Yi))v ^ 5(ri - Ti)5{r2 - r3)C(ri - rg), where 

^(r-r')-^(^)\ (4) 

g'^ \|r — r 1/ 

with g = uD the dimensionless conductance of the normal two-dimensional system. 

Indeed, a similar result was obtained by Spivak and Zhou |1^ for the perpendicular 

field geometry: the difference between the parallel and perpendicular field orientations 

is manifest only in the short-ranged Cooperon contributions which dress the vertices, 

giving rise to the envelope functions of width ~ L^ that we denoted by 6{r) — 

the long-range scaling is inherited from the field-insensitive diffuson content. While, 

for coordinates r^ separated by more than L^, matrix elements 6Il(){ri,r2) remain 

statistically uncorrelated from those of 5no(r3,r4), it may be confirmed that the 

correlation function {6Ilo{r,r)6Ilo{r',r'))v (Fig. [T|d) also scales as C(r — r'). 

Now, in the general two-dimensional system, the eigenf unctions of a typical random 

Hermitian operator exhibit a continuous evolution from strong localisation below the 

band edge of the non-disordered system, to weak localisation above. Applied to the 

pairing susceptibility, such behaviour would motivate a theory of the transition based on 

phase ordering of domains large enough to survive the effects of quantum fluctuations. 
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Figure 1. Impurity diagrams providing the leading contribution to the variance 
of the matrix element 6Ilo{r,r'). In both cases, the long-ranged character of 
the correlations can be ascribed to 'diffuson' contributions which are insensitive 
to magnetic field. The influence of the magnetic field is imposed through the 
short-ranged 'Cooperon' ladders which dress the vertices. 



However, these considerations overlook the potential significance of the long-ranged 
power-law correlations of 6Ilo{r,r'). In the present case, the complex random matrix 
elements have a strength which decays only as a power-law, with an exponent f3 that 
coincides with the dimensionality D, i.e. f3 = D = 2 in the thin-film geometry. Previous 
investigations by Levitov [HI Ej have revealed that such a dependence places the system 
in a 'critical regime' where the bulk eigenstates of the susceptibility are neither fully 
extended nor exponentially localized. Rather, such states Xa{^) exhibit a multi- fractal 
structure all the way down to the 'band edge' with moments (|xQ,(r)p^)n ~ L^'^i+'^q 
characterized by a set of exponents r^^, and power-law spatial correlations. 



fij 



L\\xi{r)xi{r')\)ii 



(5) 



At the level of the saddle-point, the transition to the hulk superconducting phase may 
not, after all, be a problem of optimal fluctuations. 

To make the analysis quantitative, one can use the variance (j3)) to characterize the 
distribution of fluctuations. On scales |r — r'| 3> Lh, let us suppose that matrix elements 
are specified by a Gaussian distribution, Ppo] = exp[— ^ J dr dr' |(5no(r, r')p/C(r — r')], 
where we include only the first type of correlations discussed above (see Fig. ^a)). 
On length scales smaller than Lh-, we will assume to be cautious that the matrix 
elements simply coincide with their impurity average. (Indeed, such an assumption 
underestimates the renormalisation of the quantum critical point.) Although there 
is no reason to expect the distribution to be Gaussian, an estimate of the leading 
contribution to the higher cumulants is compatible with the Ansatz for P[no]. With 
this distribution, the rearrangement of the DoS implied by the long-ranged fluctuations 
of Ho can be estimated within the self-consistent Born approximation (SCBA). Defining 
the impurity averaged Green function, Qq = ((7+ — n/z/)^^)n, where 7+ = 7 + iO, the 
SCBA translates to the condition 

"^/^« dp 1 



^-^o(r,r) 



(27r)2 7+-7, + L2,(pV2 + a^o) 



(6) 
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where (taking udTh ^ 1) 7^ = Ihi^ujdTh), and a = n/ q^^^qThY characterizes the 
strength of fluctuations, where k is some numerical constant. Setting a = 0, the 
self-consistent equation recovers the constant DoS of the ensemble averaged operator. 
In this approximation, as expected from the AG condition, the transition to the 
superconducting phase takes place into a state with a spatially homogeneous order 
parameter when (Abcs ^)~^ — 7c = ln(2/Aor]^) = 0. Reinstating mesoscopic fluctuations 
(through a), one finds the renormahzed edge 7c — 7c+Q;, valid in the limit (7c— 7c)/7c ^ 
1, from which one can infer the shift of the critical field. 



r'k rj, \' ■ %2 



An expansion of the DoS in the vicinity of 7c obtains p(7) = (tr(5(7 — n/z/))n — 
(3/7r^a)^/^A/7c — 7. Note that inclusion of the ((5no(r, r)5no(r', r'))y correlations leads 
only to a change in k. 

Therefore, leaving aside the potential for a small exponential band of Lifshitz tail 
states of IIo, at the level of the linearised saddle-point approximation, the onset of 
superconductivity at T = takes places at a value of I/t]^ renormahzed from that 
predicted by AG theory. The bulk transition is to a glass-like phase in which the complex 
order parameter inherits the multi-fractal structure implied by the critical theory, 

(|A(r,t)nA(r',t')|VA-(|A(r,t)p)^,^~/(r-r'). (8) 

The integrity of the saddle-point approximation depends sensitively on the impact 
of dynamical and spatial fluctuations of the order parameter. Their significance can 
be inferred from the behaviour of the impurity averaged susceptibility. The gradient 
expansion, (n^,„(r, r'))y ~ vdiy — y')[^c + L^d^, /2 — \uJm\TH/'^'T^]i shows the system to be 
dissipative (both outside and within the gapless ordered phase) with a rate proportional 
to 1/th implying a dynamical exponent z = 2. In the two-dimensional geometry, this 
places the system at its upper critical dimension jH]. In the presence of mesoscopic 
fiuctuations, the phase space for low-energy fiuctuations is only diminished justifying 
the saddle-point analysis adopted in this work. 

At temperatures T 7^ 0, the thermodynamic properties of the system demand two 
further considerations: As well as the obvious significance of thermal fiuctuations of 
the order parameter, the mechanisms of quantum interference, on which the long-range 
correlations of the pairing susceptibility rely, become gradually extinguished. Thermal 
smearing limits the long-range correlations of (5no(r,r') to length scales |r — r'| smaller 
than the thermal length Lt = y^D/T, i.e. ^(r - r') ~ {u'^/g'^){^o/\r - r'|)^e-l'-'''l/-^^. 
Providing L^ remains greatly in excess of L^, such dependence motivates a coarse- 
grained description: Over the range 1/ffj < 1/th < ^/th, each domain of size Lt has 
condensed into a glass-like superconducting droplet, weakly connected to neighboring 
domains by the residual fiuctuations of SYIq. On length scales in excess of Lt, one 
therefore expects to recover a "gauge glass" picture analogous to that envisaged by 
Spivak and Zhou [2] with a low-temperature behaviour that can be inferred from the 
analysis of Galitski and Larkin [Hj. 
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